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1 INTRODUCTION 

The TRACE observations [1-3] reveal that the solar corona is comprised of lots of thin 
loops that are intrinsically dynamic, and that continually evolve. These very thin strings, 
the observations indicate, are heated for a few to tens of minutes, after which the heating 
ceases, or at least changes significantly in magnitude [1] . In this paper we examine a class 
of mechanisms, which, through the viscous-dissipation of the plasma kinetic energy, pro- 
vide the primary and basic heating of the coronal structures during their very formation. 
The basic input of the theory is the reasonable assumption that the coronal structures 
are created from the evolution and re-organization of a relatively cold plasma flow [1-16] 
emerging from the sub-coronal region (between the solar surface and the visible corona) 
and interacting with the ambient magnetic field anchored inside the solar surface. During 
the process of trapping and accumulation, a part of the kinetic energy of the flow is con- 
verted to heat by viscous dissipation and the coronal structure is born hot and bright. For 
this to happen, we must find alternative fast and efficient heating mechanisms because, 
for the conditions prevalent in the coronal structures, the standard viscous dissipation is 
neither efficient nor fast. The rates of viscous dissipation can be considerably increased 
by processes which either enhance the local viscosity coefficient, or induce short scale 
structures in the velocity field. At present we do not know of any convincing mechanism 
for the former possibility. This paper, therefore, is limited to an examination of processes 
of the latter kind. We find that as long as the flow-velocity fleld is treated as an essential 
and integral part of the plasma dynamics, fast and desirable viscous dissipation does, 
indeed, result. Consequently, during its very formation, the coronal structure can become 
hot and bright. 

Of the several possible mechanisms by which the flow kinetic energy may be converted 
into heat we emphasize the following two: The flrst is the ability of supersonic flows to 
create nonlinear perturbations which steepen to produce short scale structures which can 
dissipate by ordinary viscosity. The second stems from the recently established property of 
the magnetofluid equilibria for extreme sub-Alfvenic flows (most of the observed coronal 
flows fall in this category) - such flows can have a substantial, fastly varying (spatially) 
velocity fleld component even when the magnetic fleld is mostly smooth. Viscous damp- 



ing associated with this varying component could be a major part of the primary heating 
needed to create and maintain the bright Corona. From a general framework describing 
a plasma with flows, we have been able to "derive" several of the essential characteristics 
of the coronal structures. Theoretical basis for both these mechanisms will be discussed. 
Our simulation (for which we developed a dissipative two-fluid code), however, concen- 
trates only on the flrst mechanism, and preliminary results reproduce many of the salient 
observational features. There is clear cut evidence of nonlinearly steepened velocity fields 
which effectively dissipate and heat the coronal structure right through the process of 
formation. The numerical investigation of the second mechanism, which will require a 
much higher spatial resolution, will be undertaken soon. 

Naturally all these processes require the existence of particle fiows with reasonable 
amounts of kinetic energy. There are several recent publications [HI] cataloguing enough 
observational evidence for such flows in the regions between the sun and the corona to 
warrant a serious investigation in this direction. It must be admitted that we still have 
little understanding of the nature of the processes by which the relatively cool material 
(no hotter than about 20000K) moves upward from low altitudes (as low as a few thou- 
sand kilometers) to the outer atmosphere. For this paper, we shall simply exploit the 
observation that the flows exist, and work out their consequences. We believe that the 
flows might prove to be a crucial element in solving the riddle of coronal heating. 

The model for the solar atmosphere that we propose and investigate is obtained by 
injecting an essential new feature into several extant notions — the plasma flows are 
allowed to play their appropriate role in determining the evolution and the equilibrium 
properties of the structures under investigation. We reiterate that the distinguishing 
ingredient of our model is the assumption (observationally suggested) that relatively cold 
particles spanning an entire range of velocity spectrum — slow as well as fast, continually 
flow from the sub-coronal to the coronal regions. It is the interaction of these cold primary 
flows with the solar magnetic flelds, and the strong coupling between the fluid and the 
magnetic aspects of the plasma that will deflne the characteristics of a typical coronal 
structure (including Coronal Holes). In this paper we limit ourselves to the formation and 
primary heating aspects; we do not deal with instabilities, their nonlinear effects, flaring 



etc. These are the problems that we will confront at the next stage of the development 
of the model. 

In Sec. 2, we desribe in relative detail our basic model for the upper solar atmosphere, 
a time-dependent, two-fluid system of currents and flows. The flows are treated at 
par with other determining dynamical quantities, the currents and the solar magnetic 
flelds. Section 3 is devoted to the derivation of the characteristics of typical coronal 
structures from the basic model. Following a general discussion, we numerically simulate 
the evolution of a cold plasma flow as it interacts with the solar magnetic fleld and gravity 
in Sec. 3.1. We follow the fate of an initial cold supersonic flow as the particles get trapped 
by the magnetic fleld. By the time a sizeable density is built up we also flnd a considerable 
rise in temperature. In a very short time the velocity fleld developes a shocklike structure 
which dissipates with ordinary viscosity to convert the flow kinetic energy to heat. In 
Sec. 3.2 we take a different approach, and describe elements of the recently investigated 
magneto-fluid theory (see Mahajan and Yoshida, 1998, 2000) which allows the existence 
of equilibrium solutions missing in the flowless MHD. We flnd that a short-scale velocity 
component is predicted to be an essential aspect of a class of magnetofluid states in terms 
of which a typical coronal structure could be modelled. The magnetofluid states are the 
equilibrium states created by the strong interaction of the magnetic and the fluid character 
of a plasma, and are derived from the normal two-fluid equations when the velocity fleld 
is treated at par with the magnetic fleld. In a somewhat detailed discussion, we argue for 
the relevance of these states for the solar corona. These states could be seen as a set of 
quasi-equilibria evolving to an eventual hot coronal structure; the dissipation of the small 
scale velocity component provides the necessary source of heating. Since the numerical 
simulation of these states requires a much flner resolution than we have in our code, their 
time dependent simulation is deferred to a future work. 



2 THE BASIC MODEL — GENERAL EQUATIONS 
FOR THE QUIESCENT SOLAR ATMOSPHERE 

In this section we will develop a general theoretical framework from which the typical 
solar coronal structure will be "derived." In our model, the plasma flows from the Sun's 
surface provide the basic source of matter and energy for the myriad of coronal structures 
(including Coronal Holes). Although the magnetic field is, naturally, the primary culprit 
behind the structural diversity of the corona, the flows (and their interactions with the 
magnetic field) are expected to add substantially to that richness. 

The primary objective of this paper is to investigate how these flows are trapped and 
heated in the closed magnetic field regions, and create one of the typical shining coronal 
elements. We shall, however, make a small digression to suggest a possible fate of the 
fast flows making their way through the regions where the magnetic field is weak, or 
has open field lines. The faster particles could readily escape the solar atmosphere in 
the open field-line regions. They could also do so by punching temporary channels in 
the neighboring closed field-line structures. The flows escaping through these existing 
or "created" coronal holes (the coronal holes (CH) are highly dynamical structures with 
open and "nearly open" magnetic field regions, see e.g. [17]) may eventually appear as 
the fast solar wind. 

In the closed field-line, the magnetic fields will trap the flows, and the trapping will 
lead to an accumulation of particles and energy creating the coronal elements with high 
temperature and density. We shall not consider the solar activity processes, since the 
activity regions (AR) and flares, though an additional source of particles and energy, 
cannot account for the continuous supply needed to maintain the corona. Moreover, in 
the theory we suggest, the flare is understood to be a secondary event and not the primary 
source for the creation of the hot corona. 

To describe the entire atmosphere of the quiescent, non-flaring Sun we use the two- 
fluid equations where we keep the flow vorticity and viscosity effects (Hall MHD). The 
general equations will apply in both the open and the closed field regions. The difference 
between various sub-units of the atmosphere will come from the initial, and the boundary 



conditions. 

Let V denotes the flow velocity field of the plasma in a region where the primary solar 
magnetic field is Bg . It is, of course, understood that the processes which generate the 
primary flows and the primary solar magnetic fields are independent (say at t = time). 
The total current j = j/ + js (here j/ is the self-current that generates the magnetic field 
Bf and j^ is the source of the solar field B^) is related to the total (that can be observed) 
magnetic field B = B^ + Bj by Ampere's law: 

j = ^VxB. (1) 

Notice that in the framework we are developing (assumption of the existence of primary 
flows), the boundaries between the photosphere, the chromosphere and the corona become 
rather artificial; the different regions of each coronal structure are distinguished by just 
the parameters like the temperature and the density. In fact, these parameters should 
not show any discontinuities; they must change smoothly along the structure. At some 
distance from the Sun's surface, the plasma may become so hot and dense that it becomes 
visible (the bright, visible corona), and this altitude could be viewed as the base of the 
corona. But to study the creation and dynamics of bright coronal structures (loops, 
arches, arcades etc.) we must begin from the photosphere, and determine the plasma 
behavior in the closed field regions. 

Assuming that the primary flows provide, on a continuous basis, the entire material 
for coronal structures, the solar flow with density n will obey the Continuity equation: 

d 

— n + V ■ (nV) = 0. (2) 

ot 

We must add a word of caution: in the closed field regions, the trapped particle density 
may become too high for the confining field, resulting in instabilities of all kinds. In this 
paper we shall not deal with instabilities and their consequences; it will constitute the 
next stage of development of the model. 

Since the corona as well as the SW are known to be mostly hydrogen plasmas (with 
a small fraction of Helium, and neutrons, and an insignificant amount of highly ionized 



metallic atoms) with nearly equal electron and proton densities: He — Ui = n , we expect 
the quasineutrality condition V ■ j = to hold. 

In what follows, we shall assume that the electron and the proton flow velocities 
are different (two-fluid approximation was used e.g. in Sturrock and Hartly, (1966). 
Neglecting electron inertia, these are Vj = V, and Vg = (V — j/en), respectively. We 
assign equal temperatures to the electron and the protons for processes associated with the 
quiescent Sun. For the creation processes of a typical coronal structure, this assumption is 
quite good. For the fast SW, however, we know from recent observations (Banaszkiewicz 
et al. 1997 and references therein), that the species temperatures are found to be different: 
Tj ~ 2 ■ 10^ K and Te ~ 1 ■ 10^ K. Since the fast SW is not the principal interest of this 
paper, we shall persist with the equal temperatures assumption; the kinetic pressure p is 
given by: 

p = Pi+p,^2nT, T = TiC^ T,. (3) 

With this expression for p, and by neglecting electron inertia, the two-fluid equations are 
obtained by combining the proton and the electron equations of motion: 

^Vu + (V ■ V)Vu = 

Vk{nT) + V, (^^) - -^ViU.^ki, (4) 
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where b = eB/mjC, -mj is the proton mass, G is the gravitational constant, Mq is the 

solar mass, r is the radial distance, and Ilj^/fc is the ion viscosity tensor. For flows with 

large spatial variation, the viscous term will end up playing an important part. To obtain 

an equation for the evolution of the flow temperature T, we begin with the energy balance 

equations for a magnetized, neutral, isothermal electron-proton plasma: 

d 

where a is the species index. The fluid energy Sa (thermal energy + kinetic energy) and 
the total pressure tensor Va^^i are given by 

^Q, = Uq I -To H — ^ 1 , Va^kl = naTa^kl + ^a,kh i^) 



and 



f« = e„E + — V„ X B + m„V- ® 



c r 

is the volume force experienced by the fluids (E is the electric field). In Eq. (6), q^ is 
the heat flux density for the species a. After standard manipulations we arrive at the 
temperature evolution equation 
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where E^ is the total radiative loss, Eh is the local mechanical heating function, and z/j is 
the ion kinematic viscosity. Note that we have retained viscous dissipation in this system. 
If primary flows are ignored in the theory, various anomalous heating mechanisms need 
to be invoked, and a corresponding term Eh has to be added. The full viscosity tensor 
relevant to a magnetized plasma is rather cumbersome, and we do not display it here. 
However just to have a feel for the importance of spatial variation in viscous dissipation, 
we display its relatively simple symmetric form. It is to be clearly understood that this 
version is meant only for theoretical elucidation and not for detailed simulation. We 
notice that even for incompressible and currentless flows, heat can be generated from the 
viscous dissipation of the flow vorticity. For such a simple system, the rate of kinetic 
energy dissipation turns out to be 
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revealing that for an incompressible plasma, the greater the vorticity of the flow, the 
greater the rate of dissipation. 

Let us now introduce the following dimensionless variables: 
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where Rq is the solar radius. Note that in general z/j is a function of density and temper- 
ature: Ui = (Vi^(/iT^/127rne^). 

In terms of these variables, our equations read: 

|V + (V.V)V^ 
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This set of equations will now be studied for different types of magnetic field regions, in 
particular the regions with closed field lines. 

Before we embark on a detailed theory of the formation and heating of the corona, we 
would like to give a short list of heating mechanisms which have been invoked to deal with 
this rather fundamental and still unresolved problem of Solar physics : Alfven waves [20- 
27], Magnetic reconnection in Current sheets [28-36], and MHD Turbulence [37-39]. For 
all these schemes, the predicted temperature profiles in the coronal structures come out to 
be highly sensitive to the form of the heating mechanism [40,41]. Parker (1988) suggested 
that the solar corona could be heated by dissipation of many tangential discontinuities 
arising spontaneously in the coronal magnetic field that is stirred by random photospheric 
footpoint motions. This theory stimulated numerous searches for observational signatures 
of nanofiares. Unfortunately, all of these attempts fall short of providing a continuous 
(both in space and time) energy supply that is required to first create in a few minutes, 
and then support for longer periods the observed bright coronal structures (see e.g. [1,2]). 



Our attempt to solve this problem makes a clean break with the conventional approach. 
We do not look for the energy source within the corona but place it squarely in the primary 
flows emerging from the Sun (see the results of [1-3]). We propose (and will test) the 
hypothesis that the energy and particles associated with the primary flows, in interaction 
with the magnetic fleld, do not only create the variety of conflgurations which constitute 
the corona, but also provide the primary heating. The flows can give energy and particle 
supply to these regions on a continuous basis — we will show that the primary heating 
takes place simultaneously with the accumulation of the corona and a major aspect of 
the flow-magnetic fleld interaction, for our system, is to provide a pathway for this to 
happen. 

A mathematical modeling of the coronal structure (for its creation and primary heat- 
ing) will require the solution of Eqs. (12)-(16) with appropriate initial and boundary 
conditions. We will use a mixture of analytical and numerical methods to extract, what 
we believe, is a reasonable picture of the salient aspects of a typical coronal structure. 



3 CONSTRUCTION OF A TYPICAL CORONAL 
STRUCTURE 

Though the solar atmosphere is highly structured, it seems that most of the constituent 
elements have something common in their creation and heating. In order to construct 
a unifled theory for the entire corona, one would have to confront large variations in 
plasma density and temperature. It seems, however, that beyond the coronal base, the 
equilibrium temperature tends to be nearly constant on each one of these structures; the 
temperature of a speciflc structure increases insigniflcantly (about 20 p.c.) from its value 
at the base to its maximum reached at the top of the structure. This change is much 
less than the temperature change (about 2 orders of magnitude) that occurs between the 
solar surface and the coronal base. This observation is an outcome of the investigation of 
several authors (see, for example, [1,2,41-45]). Their results show that the bright elements 
of the corona are composed of quasi-isothermic and ultra-thin arcs (loops) of different 
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temperature and density, situated (located) close to one other. This state is, perhaps, 
brought about by the isolating influence of magnetic fields which prevent the particle and 
energy transfer between neighboring structures. 

It is safe to assume, then, that in the quasi-equilibrium state, each coronal structure 
has a nearly constant temperature, but different structures have different characteristic 
temperatures, i.e., the bright corona seen as a single entity will have considerable temper- 
ature variation. Observations tell us that the coronal temperatures are much higher than 
those of the primary flows (which we are proposing as the mother of the corona). For the 
consistency of the model, therefore, it is essential that the primary "heating" must take 
place during the process of accumulation of a given coronal entity. 

This apparent problem, in fact, can be converted to a theoretical advantage. We 
distinguish two important eras in the life of a coronal structure; a hectic period when 
it acquires particles and energy (accumulation and heating), and the relatively calmer 
period when it "shines" as a bright, high temperature object. 

In the first era, the most important issue is that of heating while particle accumulation 
(trapping) takes place in a curved magnetic field. This is, in fact, the essential new 
ingredient of the current approach. We plan to show: 

1) that the kinetic energy contained in the primary flows can be dissipated by viscosity 
to heat the plasma, and 2) that this dissipation can be large enough to produce the 
observed temperatures. 

Naturally, a time dependent treatment will be needed to describe this era. 

Any additional heating mechanisms, operative after the emergence of the coronal 
structure, will not be discussed in this paper. For an essential energy inventory of the 
quasi-equilibrium coronal structure, we also ignore the contributions of flares and other 
"activities" on the solar surface because they do not provide a continuous and sufficient 
energy supply [2]. 

The second era is that of the quasi-equilibrium of a coronal structure of given density 
and temperature - neither of which has to be strictly constant. The primary heating has 
already been performed, and in the equilibrium state, we can neglect viscosity, resistivity 
and other collisional effects in addition to neglecting the time dependence. The calcu- 
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lations in this regime will be limited to the determination of the magnetic field and the 
velocity-field structures that the collisionless magnetofiuids can generate and we will also 
examine if these structures can confine plasma pressure. 

3.1 Creation and heating of coronal structure 

In this subsection we will concentrate on numerical methods to test our basic conjecture 
that the primary solar flows are responsible for the creation and heating of a typical 
bright coronal structure. The numerical results (obtained by modeling Eqs. (12)-(16) 
with viscosity tensor relevant to magnetized plasma) are extremely preliminary, but they 
clearly indicate that the proposed mechanism has considerable promise. 

Let us first make order of magnitude estimates on the requirements that must be 
met for this scheme to be meaningful. It is well known that (see e.g. [46]) the rate of 
energy losses F from the solar corona by radiation, thermal conduction, and advection 
is approximately 5-10^ erg/cm s. For the brightest loops the rate loss could even reach 
5 ■ 10^ erg/cm s. If the conversion of the kinetic energy in the primary fiows were to 
compensate for these losses, we would require a radial energy fiux 

]^m,n,V^ K) > F, (17) 

where Vq is the initial fiow speed. For Vq ~ 300 km/s this implies an initial density in the 
range: (3 • lO'^ - 4 • 10® )cm"^ 

For slower (~ 100 km/s) velocity primary flows the starting density has to be higher 
(> 10^ cm^^). These values seem reasonable according to the latest observational data 
[1-3]. 

The normal viscous dissipation of the flow takes place on a time (using Eq. (10)): 

^visc ~ , (18) 

where L is the length of the coronal structure. For a primary flow with Tq = 3eV = 
3.5 ■ 10^ K and no = 4 ■ 10® cm^^ creating a quiet coronal structure of size L = (2 ■ 10^ — 
7 ■ 10^°) cm, the dissipation time can be estimated to be of the order of (2 ■ 10® — 10^°) s. 
The shorter the structure and hotter the flow, the faster is the rate of dissipation. This 
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estimated time is much longer than what is actually found by the latest observations by 
TRACE [1]. Mechanisms much faster than the one embodied in (18), therefore, will be 
needed for the model to work. In the absence of "anamalous viscosity", the only way 
to enhance the dissipation rates (to the observed values) is to create spatial gradients 
of the velocity field that are on a scale much much shorter than that of the structure 
length (defined by the smooth part of the magnetic field). Thus, the viability of the 
model depends wholly on the existence of mechanisms that induce short-scale velocity 
fields. Numerical simulations show that the short-scale velocity fields are, indeed, self- 
consistently generated in the two-fiuid system. 

For numerical work (to illustrate the bright coronal structure formation), we model 
the initial solar magnetic field as a 2D arcade with circular field lines in the x-z plane 
(see Fig. 1 for the contours of the vector potential, or the fiux function). The field 
attains its maximum value B^g^^{xo, ,2 = 0) at xq at the center of the arcade, and is a 
decreasing function of the height z (radial direction). The set of model equations (12-16) 
was solved in 2D flat geometry (x,z) using the 2D version of Lax-Wendroff numerical 
scheme (Richtmyer and Morton 1967) alongwith applying the Flux-Corrected-Transport 
procedure [48]. Equation (13) was replaced with its equivalent for the y-component 
of the vector potential which automatically ensures the divergence-free property of the 
magnetic field. The equation of heat conduction was treated separately by Alternate 
Direction Implicit method with iterations [49] . Transport coefficients for heat conduction 
and viscosity were taken from Braginski, 1965. A numerical mesh of 200 x 150 points 
was used for computation. 

To illustrate the formation and heating of a general coronal structure, we have modeled 
several cases with different initial and boundary conditions for cold primary fiows. The 
dynamical picture is strongly dependent on the relation of the initial fiow pressure and the 
magnetic field strength. Two limiting cases are interesting: 1) the initial magnetic field is 
weak, and the fiow significantly deforms (and in specific cases, drags) the magnetic field 
lines, 2) the initial magnetic field is strong, and the fiow leaves the field lines practically 
unchanged. 

For sub-Alfvenic fiows, we present in Figs. 2-5 the salient features of our preliminary 
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results. We have plotted (as functions of x and z) four relevant physical quantities: the 
flux function A, the density n, the temperature T, and the magnitude of the velocity 
field |V| (for specific cases, when needed, we give the radial component of velocity field 
Vz also). 

The plots correspond to two (in some cases to three) different time frames. The results 
are described under three separate headings, covering respectively, the fully uniform, the 
spatially non-uniform, and the time-dependent as well as spatially non-uniform initial 
flows. 

A. Initially uniform primary flow and an Arcade-like magnetic field struc- 
ture — Fig. 2. 

This case is highly idealized but illustrates the main aspects of the creation of the hot 
coronal structures, and of the basic heating process. 

When discussing the temporally uniform initial flows, we choose the parameters to 
satisfy the observational constraint that, over a period of some tens of minutes, the lo- 
cation of the heating must have a relatively smooth evolution [1]. The final shape and 
location of the coronal structure (of the associated B(r,t), for example) will be naturally 
defined by its material source, by the heating dynamics, and by the initial field Bo(r,t). 

For these studies, the initial flow velocity field is taken to be uniform at the surface 
and has only a radial component, Vz = 300 km/ s. Other parameters are: Maximum value 
of the magnetic field -Bmax(a;o, -2 = 0) = TG, initial density of the flow 4 ■ 10® cm~^ and the 
initial temperature 3 eV . Simulations yield the following results: 

1) The flow particles begin to accumulate at the footpoints near the solar surface 
(Fig. 2, see density at t = 750s). The accumulation goes on with time, and gradually 
the entire volume under the arcade (starting from the central short loops) is filled with 
particles (Fig. 2, density at t = 1400 s). First the shorter loops are filled, and then the 
larger ones. 

2) The heating of the particles goes hand in hand with the accumulation (Fig. 2, plots 
for density and temperature). 

3) The regions of stronger magnetic fields are denser in population (Fig. 2, plots for 
A and n). In earlier stages of the formation of a coronal structure, the regions near the 
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base (where the field is stronger) are denser and hotter than the distant regions (Fig. 2, 
t = 750 s plots for n and T ); for shorter loops, the density increases (as a function 
of height z) from the bottom of the structure, and then falls — first rapidly and later 
insignificantly; the maximum density is much greater than the initial density of the fiows. 

4) The dissipation of the fiow kinetic energy is faster in the first stage of formation 
(Fig. 2, t = 750s plot for \V\ ). The plot \V\ versus z shows steep (shock-like) gradients 
near the base. Thus the bright base is created in the very first stage in the stronger 
magnetic field regions (shorter loops). For given parameters, the initial flow is strongly 
supersonic. Thus the shocks are generated with efficient transfer of kinetic energy into 
heat. As the mean free path of ions in the plasma is of the order of (10^ — 10"^ )cm (in 
the direction parallel to the magnetic field) and the dimension of the structure is much 
greater - of the order of 10^^ cm - efficient conditions for the kinetic energy dissipation 
exist. The plots for the velocity, temperature and density reveal that with increasing z, 
and in the regions away from the arcade center, we first find an undisturbed fiow with low 
temperature, then see a transient area with high density and temperature, and finally a 
shock consistent with Hugoniot conditions. The short scale represented by the width of 
the shock-layer (determined by viscosity) is the main enhancer of viscous dissipation. 

5) For later times, the brightening process spreads over wide regions (Fig. 2, t = 1400 s 
plot for temperature). 

6) In the very first stage, the shorter loops are a bit overheated, but they cool down 
somewhat at later times when the longer loops begin to get heated (Fig. 2, plots for 
temperature). 

7) The base (T > 100 eV) of the bright region is at about 1.4-109 cm ~ O.O2i?0 (Fig. 2, 
t = 1400 s plots for n and T) from the solar surface. This number is in a very good 
agreement with the latest TRACE results [1]. Outwards from the base, the accumulated 
layer has somewhat lower, but more or less uniform, insignificantly decreasing density. In 
the accumulated layer the kinetic energy of the fiow is essentially uniform (again, decreases 
insignificantly); the dissipation has practically stopped (Fig. 2, t = 1400 s = 23mm , 
plot for \V\ versus z). The temperature is practically uniform in the longer loops and 
increases insignificantly in shorter loops (for some special conditions these conclusions 
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may be somewhat modified in specific regions of the arcade; see point 8)). Outwards 
from the hottest region of the arcade, the temperature decreases gradually and at some 
radial distance the outer boundary of the bright part is reached (Fig. 2, t = 1400 s plot 
for temperature). Thus, in a very short time a dense and bright "coronal structure" is 
created — this object survives for a time much longer than was needed for its creation. The 
simulations show that the heating process may continue during this so-called equilibrium 
stage, but at a rate much slower than the earlier primary heating. This heating seems 
just additional and supporting to the heat content of the nascent hot structure. At this 
time, however, the velocity field is already much smaller in magnitude as compared to 
the initial values; the flows in the hot coronal structure are already subsonic. This is a 
possible explanation why supersonic flows may not be seen in the hot observable coronal 
structures. 

8) When relatively dense primary flows interact with weak arcade-like magnetic fields 
{B„iax{xo,ZQ = 0) < 10 G for our initial flow with given above parameters), the field 
lines begin to deform (soon after the creation of the solar base) in the central region of 
the arcade but far from the base (see t = 1400 s plots for density and temperature in 
Fig. 2). The particle accumulation is still strong, and the dissipation, though quite fast, 
stops rather rapidly. Consequently, the temperature first reaches a maximum (up to the 
deformed field-line region this maximum is reached at the summit for each short loop) 
and later falls rapidly. Gradually one can see signs for the creation of a local gravitational 
potential well behind the shortest loops (see t = 1400 s plot for A in Fig. 2). This well 
supports a relatively dense and cold plasma in the central area of the arcade (t = 1400 s 
for n and T of Fig. 2). The density of this structure is considerably greater than that of 
the surrounding areas, and the temperature is considerably less than that of the rest of 
the accumulated regions at the same height of the arcade. 

Our preliminary simulations show that for the same parameters of the primary outflow, 
such cold and dense plasma objects (confined in the so-called potential well) will not form 
in the regions where the initial magnetic fields are stronger (i?max(a;o5 -Zo = 0) > 20 G). 

B. Spatially non uniform primary flow interacting with an arcade— like mag- 
netic field structure — Figs. 3, 4. 
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The latest observations support the idea that the coronal material is injected disconti- 
nously (in pulses or bunches, for example) from lower altitudes into the regions of interest 
(e.g. spicules, jet-like structures [6,7,12,13,1,2]). A realistic simulation, then, requires a 
study of the interaction of spatially non-uniform initial flows with arcade-like magnetic 
field structures. These "close to the actual" cases represent more vividly the dynamics of 
the hot coronal formation. 

1) When the spatially symmetric initial flow (plot for \4 at t = in Fig. 3a) inter- 
acts with the arcade (plots in Fig. 3), and the initial magnetic field is rather strong 
{Bi-aax{xo, zq = 0) = 20 G) , the primary heating is completed in a very short time 
(~ (2 — 3)min) on distances (~ 10000 km) shorter than the uniform-flow case when the 
initial magnetic fleld was weaker. This is also consistent with observations. The heating 
is very symmetric and the resulting hot structure is uniformly heated to 1.6 • lO^K. 

2) Observations reveal the existence of cool material and downflows, right within the 
hot coronal structures; they also show an imbalance in the primary heating on the two 
sides of the loops (see [12, 1]). To reproduce these characteristics, we have modeled the 
coronal structure formation process using an asymmetric, spatially non-uniform initial 
flow interacting with a strong magnetic fleld (see Fig. 4). 

For both of the discussed cases, the downflows can be clearly seen for the velocity fleld 
component Vz- In Fig. 4, the downflow is created simply by changing the initial character 
of the flow (initially we had only the right pulse from the velocity fleld distribution given 
in Fig. 3a), while in Fig. 3a (plot at t = 297 s), the downflows are the result of more 
complicated events (see explanation below, in the next paragraph). The flnal parameters 
of the downflows are strongly dependent on the initial and boundary conditions. In the 
pictures, the imbalance in the primary heating process is also revealed. 

When two identical pulses (Fig. 3a, plot at t = s) enter in succession into our 
standard, arcade-like initial magnetic fleld, we simulate the equivalent of two colliding 
flows on the top of a structure. Shocks, though not very strong, are generated in a very 
short time (t = 30 s). Such shocks, on both sides of the arcade-center, have hot fronts 
and cold tails. Soon (t = 42s) these shocks become "visible", a hot and dense area is 
created on top of the structure where these shocks (at this moment they have become 
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stronger) collide. After the collision (and "reflection"), the entire area within the arcade 
becomes gradually hot. At some moment, a practically uniformly heated structure is 
created, and the primary heating stops. This process is accompanied by downfiows much 
slower than the primary flows; much of the primary flow kinetic energy has been converted 
to heat via shock generation (the shock and downfiow velocities differ significantly). It 
is clear that in the case of spatially assymetric initial flows, the downfiows on different 
sides of the arcade-center will have different characteristics. Due to the high pressure 
prevalent in the nascent hot structure (loop), there is no more inflow of the plasma and 
the flow deposits its energy at the base; the base becomes overheated. Later this energy 
can be again transferred upwards via thermal conduction (this mechanism can work in 
all the discussed cases), but at that moment the flow could be also changed (see initially 
time-dependent flow cases below). 

Plots for the temperature and velocity field in Figs. 3b, 4 also indicate that some cold 
particles still remain in the body of the newly created hot structure. These particles are 
perhaps from the slower aggregates (our initial flow was not uniform) which did not have 
sufficient energy to be converted to heat. 

C Time dependent non uniform initial flows interacting with arcade— like 
magnetic field structures. Fig. 5 

To simulate reality further we introduce time dependence in the initial primary flow 
velocity field. We discuss two distinct cases: 

f ) Initially, the velocity field has a pulse-like distribution with a time-period nearly 
half of the "formation time" of the quasi-equilibrium structure corresponding to the case 
with time-independent initial conditions. The results displayed in Fig. 5 show that the 
emerging coronal structure has a rather uniform distribution of temperature along the 
magnetic field, and the latter is practically undeformed during formation and heating. 
We see that when the basic heating ceases, the hot structure survives for the time of 
computation which happens to be shorter than the time necessary for losses that destroy 
the structure. 

2) The velocity field has a fast amplitude modulation near its maximum value ( for 
these simulations the maximum radial velocity was taken to be 300A;m/s). We find 
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that the dynamics of the hot coronal structure creation is quite similar to the initially 
time-independent, spatially symmetric case. Because of this, we don't give here the 
corresponding plots. We only note that for this case, the structure tends to become 
even hotter (by a factor 1.2 for the same parameters) and when quasi-equilibrium is 
established (time for this to happen is longer than for the time-independent initial flows) 
the base of the structure is hotter than the top although at an earlier time the top was 
hotter, i.e, there is a temperature oscillation with a time-period longer compared to the 
creation time of the hot structure. 

The main message of numerical simulation is that the dynamical interaction of an 
initial flow with the ambient solar magnetic field leads to a re-organization of the plasma 
such that the regions in the close vicinity of the solar surface are characterized by strongly 
varying (in space and time) density and temperature, and even faster varying velocity 
field, while the regions farther out from the bright base are nearly uniform in these 
physical parameters. This phenomenon pertains generally, and not for just a set of specific 
structures. The creation and primary heating of the coronal structures are simultaneous, 
accompanied by strong shocks. These are fast processes (few tens of minutes) taking 
place at very short radial distances from the Sun (~ 10000 km) in the strong magnetic 
field regions with significant curvature. The final characteristics of the created coronal 
structures are defined by the boundary conditions for the coupled primary fiow-solar 
magnetic field system. The stronger the magnetic field, the faster is the process of creation 
of the hot coronal structure with its base nearer the solar surface. To investigate the near 
surface region one must use general time-dependent 3D equations. Quasi-stationary 
(equilibrium) equations, on the other hand, will suffice to describe the hot and bright 
layers — the already existing visible coronal structures. 

3.2 Construction of quasi— equilibrium coronal structure 

The familiar MHD theory (single-fiuid) is a reduced case of the more general two-fluid 
theory discussed in this paper. Constrained minimizaion of the magnetic energy in MHD 
leads to force-free static equilibrium configurations [50,51]. The range of two-fluid relaxed 
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states, however, is considerably larger because the velocity field, now, begins to play an 
independent fundamental role. The presence of the velocity field not only leads to new 
pressure confining states [52,53], but also to the possibility of heating the equilibrium 
structures by the dissipation of kinetic energy. The latter feature is highly desirable if 
these equilibria were to be somehow related to the bright coronal structures. 

We begin investigating the two-fiuid states by first studying the simplest, almost 
analytically tractable, equilibria. This happens when the pressure term in the equation 
of motion (12) becomes a full gradient, i.e, whenever an equation of state relating the 
pressure and density can be invoked. For our present purpose, we limit ourselves to the 
constant temperature states allowing n~^Vp -^ 2 TV In n . 

Normalizing n to some constant coronal base density no (reminding the reader that 
no is different for different structures!), and using our other standard normalizations 
(Ajo = c/uio is defined with no), our system of equations reduces to: 

-V X b X b + V I — - /3o In n - — 1 + V X (V X V) = 0, (19) 

n \ r 2 J 

V X (v - —V X b") X b = 0, (20) 

V • (nV) = 0, (21) 

where r^o, tto, Po are defined with no, Tq, Bq. This is a complete system of seven 
equations in seven variables. 

Following Mahajan and Yoshida (1998) and [54], we seek equilibrium solutions of the 
simplest kind. Straightforward algebra leads us to the following system of linear equations: 

b + aoV x\ = dn\ (22) 

and 

b = an [v-— Vxbl, (23) 

n J 

where a and d are dimensionless constants related to the two invariants: the magnetic 
helicity /(A • B) d^x and the generalized helicity /(A + V) ■ (B + V x Y)d^x (or 
/(V-B + A-VxV + V-VxV) d^x ) of the system. We will discuss a and d later. The 
equilibrium solutions (22), (23) encapsulate the simple physics: 1) the electrons follow 
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the field lines, 2) while the ions, due to their inertia, follow the field lines modified by the 
fluid vorticity. These equations, when substituted in (19), (20), lead to 

v(^-Po\nn-^)=0, (24) 

giving the Bernoulli condition which will determine the density of the structure in terms of 
the flow kinetic energy, and solar gravity. Equations (22) and (23) are readily manipulated 
to yield 

^'^-V xVxV + aoVx f--rfnW+(l--)v = 0. (25) 



n \a / \ a 

which must be solved with (24) for n and V; the magnetic field can, then, be determined 

from (22). 

Equation (24) is solved to obtain 



n = exp 






(26) 



where g{r) = Vco/r. This relation is rather interesting; it tells us that the variation in 
density can be quite large for a low /Sq plasma (coronal plasmas tend to be low (3q; the 
latter is in the range 0.004 — 0.05) if the gravity and the flow kinetic energy vary on length 
scales comparable to the extent of the coronal structure. In this system of equations, as 
we mentioned above, the temperature (which defines /3o) has to be fixed by initial and 
boundary conditions at the base of the structure. Substituting (26) into (25) will yield a 
single equation for velocity which is quite nontrivially nonlinear. Numerical solutions of 
the equations are tedious but straightforward. 

For analytical progress, essential to revealing the nature of the self-consistent fields 
and flows, we will now make the additional simplifying assumption of constant density. 
This is a rather drastic step (in numerical work, we take the density to be a proper 
dynamical variable) but it can help us a great deal in unraveling the underlying physics. 
There are two entirely different situations where this assumption may be justified: 

1) the primary heating of corona has already been performed, i.e., a substantial part of 
flow initial kinetic energy has been converted to heat. The rest of the kinetic energy, i.e., 
the kinetic energy of the equilibrium coronal structure is not expected to change much 
within the span of a given structure. Note that the ratio of velocity components will 
have a large spatial variation, but the variation in V'^ is expected to be small. It is also 
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easy to estimate that within a typical structure, gravity varies quite insignificantly. There 
will be exceptional cases like the neighborhood of the Coronal holes and the streamer 
belts, where significant heating could still be going on, and the temperature and density 
variations could not be ignored. Such regions are extremely hard to model; 

2) if the rates of kinetic energy dissipation are not very large, we can imagine the 
plasma to be going through a series of quasi-equilibria before it settles into a particular 
coronal structure. At each stage we need the velocity fields in order to know if an appro- 
priate amount of heating can take place. The density variation, though a factor, is not 
crucial in an approximate estimation of the desired quantities. 

The constant density assumption n = 1 will be used only in Eq. (25) to solve for the 
velocity field (or the b field which will now obey the same equation). These solutions, 
when substituted in Eq. (26), would determine the density profile (slowly varying) of a 
given structure. 

In the rest of this sub-section we will present several classes of the solutions of the 
following linear equation: 

a^ V X V X Q + «o (- - d]\/ X Q + ll - -] Q = 0, (27) 

where Q is either V or b. To make contact with existing literature, we would use b as our 
basic field to be determined by Eq. (25); the velocity field V will follow from Eqs. (22) 
and (23), which for n = 1, become 

b + aoV X V = rfV (22') 

and 

b = a [V - ao V X b] . (23') 

It is worth remarking that in order to derive the preceding set of equations, all we 
need is the constant density assumption; the temperature can have gradients and, these 
are determined from the Bernoulli condition (20) with (3o{T) replacing /5oln n. 
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3.3 Analysis of the Curl Curl Equation, Typical Coronal Equi- 
libria 

The Double Curl equation (27) was derived only recently [53] (Mahajan and Yoshida 
1998); its potential, is still, largely unexplored(see [53], [54]). The extra double curl (the 
very first) term distinguishes it from the standard force-free equation [55,50,56] (Woltjer 
1958; Taylor 1974, 1986; Priest 1994 and references therein) used in the solar context. 
Since a and d are constants, Eq. (25), without the double curl term, reproduces what has 
been called the "relaxed state" [50,56]. We will see that this term contains quantitative 
as well as qualitative new physics. 

In an ideal magnetofiuid, the parameters a and d are fixed by the initial conditions; 
these are the measures of the constants of motion, the magnetic helicity, and the fiuid 
plus cross helicity or some linear combination thereof [53,52,57,38]. In our calculations, 
a and d will be considered as given quantities. The existence of two, rather than one (as 
in the standard relaxed equilibria) parameter in this theory is an indication that we may 
have, already, found an extra clue to answer the extremely important question: why do 
the coronal structures have a variety of length scales, and what are the determinants of 
these scales? 

We also have the parameter a^, the ratio of the ion skin depth to the solar radius. For 
typical densities of interest (~ (10^ — 10^)cm~^), its value ranges from (~ 10"'^ — 10^^); 
a very small number, indeed. Let us also remind ourselves that the |V| is normalized to 
the inverse solar radius. Thus | V| of order unity will imply a structure whose extension 
is of the order of a solar radius. To make further discussion a little more concrete, let us 
suppose that we are interested in investigating a structure that has a span ei?©, where e 
is a number much less than unity. For a structure of order 1000 km, e ~ 10^^. The ratio 
of the orders of various terms in Eq. (25) are (|V| ~ L^^) 

4: 20 fi_rfV h _ 



(2^ 
(1) (2) (3) 

Of the possible principal balances, the following two are representative: 
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(a) The last two terms are of the same order, and the first <^ them. Then 

e ~ «o / ,, ■ 29 

1 — a/ a 



For our desired structure to exist (ao ~ 10 ^ for no ~ 10^ cm ^), we must have 

10^ (30) 



^/a-d ^ . 



1 — d/a 

which is possible if d/a tends to be extremely close to unity. For the first term to be 
negligible, we would further need 

— <--c^^e>T7 ^, (31) 

e a 1/a — d 

which is easy to satisfy as long as neither of a ~ (i is close to unity. This is, in fact, the 
standard relaxed state, where the flows are not supposed to play an important part for 
the basic structure. For extreme sub-Alfvenic flows, both a and d are large and very close 
to one another. Is the new term, then, just as unimportant as it appears to be? The 
answer is no; the new term, in fact, introduces a qualitatively new phenomenon: Since 
V X (V X b) is second order in |V|, it constitutes a singular perturbation of the system; 
its effect on the standard root (2) ~ (3) ^ (1) will be small, but it introduces a new root 
for which the |V| must be large corresponding to a much shorter length scale (large |V|). 
For a and d so chosen to generate a 1000 km structure for the normal root, a possible 
solution would be d/a ~ 1 + 10~^, d c:^ a = —10 , then the value for |V| for the new root 
will be (the balance will be from the flrst two terms) 



|Vr^ ~ 10^ cm, 

that is, an equilibrium root with variation on the scale of 100 cm will be automatically 
introduced by the flows. The crucial lesson is that even if the flows are relatively weak 
(a ~ (i ~ 10), the departure from Vx B = aB, brought about by the double curl term can 
be essential because it introduces a totally different and small scale solution. The small 
scale solution could be of fundamental importance in understanding the effects of viscosity 
on the dynamics of these structures; the dissipation of these short scale structures may 
be the source of primary plasma heating. 
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We do understand that to properly explain the parallel (to the field-line) motion one 
must use kinetic theory since the mean free path along B lines can become of the order 
of (10^ — 10^)cm for the hot plasma (100 eV). But since the dissipation acts on the 
perpendicular energy of the fiow, we expect the two-fluid theory to give qualitatively 
(and even quantitatively) correct results. 

We would like to remind the reader that by manipulating the force free state V x B = 
a(x)B, Parker has built a mechanism for creating discontinuities (short scales) (Parker 
1972, 1988, 1994). It is important to note that short length scales are automatically there 
if plasma flows are properly treated. 

(b) The other representative balance arises when we have a complete departure from 
the one-parameter, conventional relaxed state. In this case, all three terms are of the 
same order. In the language of the previous section, this balance would demand 

1 l/a-d 

e ~ «o —, -1 ~ tto -TT- (32) 

i/a — a 1 — a/ a 



which translates as: 



and 



--d] r^l-- (33) 

a J a 



(i ~ ao -• (34) 

a e 



For our example of a 1000 km structure, a^ ■ 1/e ~ 10~^, both a and d not only have to 
be awfully close to one another, they have to be awfully close to unity. To enact such a 
scenario, we would need the flows to be almost perfectly Alfvenic. However, let us think 
of structures which are on the km or 10 km size. In that case a^ ■ 1/e ~ 10^^ or 10~^, 
and then the requirements will become less stringent, although the flows needed are again 
Alfvenic. At a density of (1 — 4) ■ lO^cm"^, and a speed ~ (200 — 300) km/s, the flow 
becomes Alfvenic for i?o ~ (1 — 3) G. It is possible that the conditions required for such 
flows may pertain only in the weak magnetic fleld regions. 

Following are the obvious characteristics of this class of flows: 

(1) Alfvenic flows are capable of creating entirely new kinds of structures, which are 
quite different from the ones that we normally deal with. Notice that here we use the term 
flow to denote not the primary emanations but the plasmas that constitute the existing 
coronal structures, or the structures in the making. 
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(2) Though they also have two length scales, these length scales are quite comparable 
to one another: This is very different from the extreme sub-Alfvenic flows where the 
spatial length-scales are very disparate. 

(3) In the Alfvenic flows, the two length scales can become complex conjugate, i.e., 
which will give rise to fundamentally different structures in b and V. 

Defining p = (1/a — d) and g = (1 — d/a), Eq. (27) can be factorized as 

(aoV X -A)(aoV x -fi)h = (35) 

where A(A+) and /i(A_) are the solutions of the quadratic equation 



p p"^ 
ao\± = --±\J—-q. (36) 

If Ga is the solution of the equation 

V X G(A) = AG(A), (37) 

then it is straightforward to see that 

h = axG{X) + a^G{fi), (38) 

where ax and a^ are constants, is the general solution of the double curl equation. Using 
Eqs. (23'), (37), and (38), we find for the velocity field 

V = - + aoV X b = ( - + aoX] axG(\) + ( - + ao/i ) anG(fi). (39) 

a \a / \a / 

Thus a complete solution of the double curl equation is known if we know the solution 
of Eq. (37). This equation, also known as the 'relaxed-state', or the constant A Beltrami 
equation, has been thoroughly investigated in literature (in the context of solar astro- 
physics see for example Parker (1994); Priest (1994)). We shall, however, go ahead and 
construct a class of solutions for our current interest. The most important issue is to be 
able to apply boundary conditions in a meaningful manner. 

We shall limit ourselves to constructing only two-dimensional solutions. For the Carte- 
sian two-dimensional case {z representing the radial coordinate and x representing the 
direction tangential to the surface, d/dy = 0) we shall deal with sub-Alfvenic solutions 
only. This is being done for two reasons: 1) The flows in a majority of coronal structures 
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are likely to be sub-Alfvenic, and 2) this will mark a kind of continuity with the literature. 
The treatment of Alfvenic flows will be left for a future publication. 

We recall from earlier discussion that extreme sub-Alfvenic flows are characterized by 
a ~ c? ^ 1. In this hmit, the slow scale A ~ (li — a)/aoda, and the fast scale /x = d/ao, 
and the velocity field becomes 

y = -axGx + daf,G{fi) (40) 

revealing that, while, the slowly varying component of velocity is smaller by a factor 
(a^^ ^ d^^) as compared to the similar part of the magnetic field, the fast varying 
component is a factor of d larger than the fast varying component of the magnetic field! 
In a magnetofiuid equilibrium, the magnetic field may be rather smooth with a small 
jittery (in space) component, but the concomitant velocity field ends up having a greatly 
enhanced jittery component for extreme sub-Alfvenic fiows (Alfven speed is defined w.r. 
to the magnitude of the magnetic field, which is primarily smooth, and for consistency we 
will insure that even the jittery part of the velocity field remains quite sub-Alfvenic). We 
shall come back to elaborate this point after deriving expressions for the magnetic fields. 
Equation (37) can also be written as 

V^G{X) + X^G{X) = 0, (41) 

and solving for one component of G(A) determines all other components up to an inte- 
gration. For the boundary value problem, we will be interested in explicitly solving for 
the z (radial) component. 

The simplest illustrative problem we solve is the boundary value problem in which we 
specify the radial magnetic field bz{x,z = 0) = f{x), and the radial component of the 
velocity field Vz{x,z = 0) = vog{x), where vq (~ d^^ <^ 1) is explicitly introduced to 
show that the fiow is quite sub-Alfvenic. A formal solution of (^^(A) = Q\) 

^,^,^g,.0 ,42) 

ox^ oz^ ag 

may be written as 

Q^= r dke-^^' Cke'''''+ f '''' dk cosqxz Ake''''' + c.c. (43) 

J\/ao Jo 
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where kx = (/c^ — X^/alY^"^, q\ = (A^/q;q — A;^)^/^, and Ck and Ak are the expansion 
coefficients. The equivalent quantities for Q^ are k^, g^, Dk-, and ii^fc. The boundary 
conditions at 2; = yield (we absorb an overall constant in the magnitude of 6^, and 
a^/ax is absorbed in Dk and Ek): 

f{x) = Qx{z = 0) + Q^{z = 0), (44a) 

vo g{x) = - Qx{z = 0) + d Q^{z = 0). (44b) 

Taking Fourier transform (in x) of Eq. (44), we find, after some manipulation, that {vq ~ 
d-\ \m\ c:^ \g{k)\) 

Ck ^ m, (45) 



-.- ^2 +^9{k)^d-'f{k), (46) 

and functionally (in their own domain of validity) Ck = Ak and Dk = Ek- With the 
expansion coefficients evaluated in terms of the known functions (their Fourier transforms, 
in fact), we have completed the solution for bz, Vz and hence of all other field components. 
The most remarkable result of this calculation can be arrived at even without a nu- 
merical evaluation of the integrals. Although /(A;) and ^(/c) are functions, we would 
assume that they are of the same order |/(/c)| = |^(A;)|. Then for an extreme sub-Alfvenic 
fiow (|V| ~ d~^ ~ 0.1, for example), the fastly varying part of hz[Q^ is negligible 
(~ d^'^ = 0.01) compared to the smooth part (Qx)- However, for these very parameters, 
the ratio 

K(A) \dDk\ \Ck/d\ ' ^ ' 

the velocity field is equally divided between the slow and the fast scales. We believe 
that this realization may prove to be of extreme importance to Coronal physics. Viscous 
damping of this substantially large as well as fastly varying fiow component may provide 
the bulk of primary heating needed to create and maintain the bright, visible Corona. 

The preceding analysis warns us that neglecting viscous terms in the equation of 
motion may not be a good approximation until a large part of the kinetic energy has been 
dissipated. It also appears that the solution of the basic heating problem may have to 
be sought in the pre-formation rather than the post-formation era. Our time dependent 
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numerical simulation to study the formation of coronal structures was strongly guided by 
these considerations. 

It is evident that for extreme sub-Alfvenic flows, the magnetic field, unlike the velocity 
field, is primarily smooth. But for strong flows, the magnetic fields may also develop 
a substantial fastly varying component. In that case the resistive dissipation can also 
become a factor to deal with. We shall not deal with this problem in this paper. 

Depending upon the choice of f{x) (from which f{k) follows) we can construct loops, 
arcades and other structures seen in the corona. 



4 Conclusions and discussions 

In this paper we have investigated the conjecture that the structures which comprise the 
solar corona (for the quiescent Sun) owe their origin to particle (plasma) flows which enter 
the "coronal regions" from lower altitudes. These primary emanations (whose eventual 
source is likely to be the sun itself) provide, on a continuous basis, much of the required 
material and energy which constitutes the corona. From a general framework describing 
a plasma with flows, we have been able to "derive" several of the essential characteristics 
of the typical coronal structures. 

The principal distinguishing component of the investigated model is the full treatment 
accorded to the velocity fields associated with the directed plasma motion. It is the 
interaction of the fluid and the magnetic aspects of the plasma that ends up creating so 
much diversity in the solar atmosphere. 

This study has led to the following preliminary results: 

1. By using different sets of boundary conditions, it is possible to construct various 
kind of 2D loop and arcade configurations. 

2. In the closed magnetic field regions of the solar atmosphere, the primary flows can 
accumulate, in periods of a few minutes, sufficient material to build a coronal structure. 
The ability of the supersonic flows to generate shocks, and the viscous dissipation of these 
shocks can provide an efficient and sufficient source for the primary plasma heating which 
may take place simultaneously with the accumulation. The stronger the spatial gradients 
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of the flow, tlie greater is the rate of dissipation of the kinetic energy into heat. The hot 
base of the structures is reached at typical distances of a ~ 10000 km from the origin of 
simulation. 

3. A theoretical study of the magnetofluid equilibria reveal that for extreme sub- 
Alfvenic flows (most of the created corona flows) the velocity field can have a substantial, 
fastly varying (spatially) component even when the magnetic field may be mostly smooth. 
Viscous damping associated with this fast component could be a major part of the pri- 
mary heating needed to create and maintain the bright, visible coronal structure. The 
far-reaching message of the equilibrium analysis is that neglecting viscous terms in the 
equation of motion may not be a good approximation until a large part of the kinetic 
energy in the primary flow has been dissipated. 

4. The qualitative statements on plasma heating, made in points 1 and 2, were tested 
by a numerical solution of the time-dependent two-fluid system. For sub-Alfvenic pri- 
mary flows we flnd that the particle-accumulation begins in the strong magnetic fleld 
regions (near the solar surface), and soon spans the entire volume of the closed magnetic 
field region. It is also shown that, along with accumulation, the viscous dissipation of the 
kinetic energy contained in the primary flows heats up the accumulated material to the 
observed temperatures, i.e., in the very flrst (and fast, ~ (2 — 10) min) stage of accumu- 
lation, much of the flow kinetic energy is converted to heat. This happens within a very 
short distance (transition region) of the solar surface ~ O.O3-R0. In the transition region, 
the flow velocity has very steep gradients. Outside the transition layer the dissipation is 
insigniflcant, and in a very short time a nearly uniform (with insigniflcantly decreasing 
density and temperature on the radial distance), hot and bright quasi-equilibrium coronal 
structure is created. In this newborn structure, one flnds rather weak flows. One also 
finds downfiows with their parameters determined by the initial and boundary conditions. 

The transition region from the solar surface to this equilibrium coronal structure is 
also characterized by strongly varying (both radial and across) temperature and density. 
Depending on the initial magnetic field , the base of the hot region (of the bright part) 
of a given structure acquires its appropriate density and temperature. 

5) The details of the ensuing dynamics are strongly dependent on the relative values of 
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the pressure of the initial flow, and of the ambient solar magnetic field in the region. Two 
limiting cases were studied with the expected results: 1) The flow entering a relatively 
weak initial magnetic field strongly deforms (and in specific cases drags) the magnetic 
field lines, and 2) the flow interacting with a relatively strong magnetic field leaves it 
virtually unchanged. 

We end this paper with several qualifying remarks: 

1) This study, in particular the numerical work, is preliminary. We hope to be able to 
extend the numerical work to make it considerably more quantitative, and to cover a much 
greater variety of the initial and boundary conditions to simulate the immense coronal 
diversity. Then a thorough comparison with observations can be undertaken. To show 
the dissipation of small scale velociy component just like the dissipation of shock-like 
structures is postponed for future since it requires much higher resolution. 

2) This paper is limited to the problem of the origin, the creation and the primary 
heating of the coronal structures. The processes which may go on in the already existing 
bright equilibrium corona (secondary or supporting heating, instabilities, reconnection) 
etc., for example, are not considered. Because of this lack of overlap between our model 
and the conventional coronal heating models, we do not find it meaningful to compare 
our work with any in the vast literature on this subject. Led by observations alone, we 
have constructed and investigated the present model. 

3) We do not know much about the primary solar emanations on which this entire 
study is based. The merit of this study, however, is that as long as they are present (see 
e.g. [1-3]), the details about their origin are not crucial. 

4) We are just beginning to derive the consequences of according a co-primacy (with 
the magnetic field) to the flows in determining overall plasma dynamics. The addition 
of the velocity fields (even when they are small) brings in essential new physics, and will 
surely help us greatly in understanding the richness of the plasma behavior found in the 
solar atmosphere. 
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Figure Captions 

Fig. 1 

Contour plots for the vector potential A (flux function) in the x — z plane for a 
typical arcade-like solar magnetic field (initial distribution). The field has a maximum 

Bmax{Xo = 0, Zo = 0) = 7G . 

Fig. 2 

Hot coronal structure formation by the interaction of the spatially homogenuous pri- 
mary flows with 2D arcade-like structure given in Fig. 1 . The initial parametrs are: 
Vzo = 300 km/s, the temperature and density of the flow, Tq = 3 eV and no = 4 ■ 10^ cm~^ 
respectively, and the background density = 10^ cm^^. The vector potential A, the flow 
density n (normalized to no), the flow temperature T (in eV) and the magnitude of the 
flow velocity \V\ (in cm/s) are plotted for t = 750s and t = 1400s . The base of the 
hot structure is created at a radial distance ~ 14000 km. The distnace scale on the plots 
is 1 = 4 ■ lO^^cm . The primary heating (and brightening) of the structure is practically 
stopped in about 23 minutes. 

Fig. 3a 

The distribution of the radial component Vz (with a maximum of 300 km/s at t = ) 
for the symmetric, spatially non-uniform velocity field . The plot scale is 1 = 5 ■ 10^ cm,. 
The process of interaction of such primary flows with the arcade-like magnetic fields 
(given in Fig. 1 with Bmax = 20 G) is accompanied by downflows much slower than the 
primary flows (plot for Vz at t = 297 s). The final parameters of downflows are strongly 
dependent on the initial and boundary conditions. 

Fig. 3b 

Hot coronal structure formation by the interaction of the initially symmetric spatially 
non-uniform primary flows (see plot for Vz{x,z) in Fig. 3a) with the 2D arcade-like 
structure given in Fig. 1 . Initial parameters are: the temperature and density of the flow. 
To = 3 eV and hq = 4-10^ cm,^^ respectively, the initial background density = 2 ■ 10^ cm,~^, 
and the fleld maximum Bmax{xo, -^o = 0) = 20 G. The plot scale is 1 = 5 ■ 10^ cm. The 
primary heating is completed in a very short time ~ (2—3 )m,in on distances (~ 10000 km,) 

38 



shorter than the uniform-flow case when magnetic field was weaker. The heating is 
symmetric and the resulting hot structure is uniformly heated to 1.6 ■ 10^ K . Much of the 
primary flow kinetic energy has been converted to heat via shock generation. 

Fig. 4 

The interaction of an initially asymmetric, spatially non- uniform primary flow (just 
the right pulse from the distribution given in Fig. 3a ) with a strong arcade-like magnetic 
field {Bmaxi^O) -2o = 0) = 20 6"). Downflows, and the imbalance in primary heating are 
revealed. 

Fig. 5 

The interaction of the time-dependent non-uniform initial flow (see plot for the time- 
distribution of Vz in this Figure; the spatial distribution of the pulse is the same as 
in Fig. 3a) with the arcade-like magnetic field structure (plot in Fig.l with B^ax = 
20 G). The emerging coronal structure has uniform distribution of temperature along the 
magnetic field (plot for T dX t = 371 s ) and the latter is practically undeformed during 
the formation and heating. 
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